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Abstract 



Power counting is applied to relativistic mean-field energy functionals to 
estimate contributions to the energy from individual terms. New estimates for 
isovector, tensor, and gradient terms in finite nuclei are shown to be consistent 
with direct, high-quality fits. The estimates establish a hierarchy of model 
parameters and identify how many parameters are well constrained by bulk 
nuclear observables. We conclude that four (possibly five) isoscalar, non- 
gradient parameters, one gradient parameter, and one isovector parameter 
are well determined by the usual bulk nuclear observables. 
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I. INTRODUCTION AND MOTIVATION 



Relativistic mean-field (RMF) descriptions of tlie nuclear many-body problem have ex- 
isted for nearly fifty years |Jl|] . An extensive body of work pj-Q has shown that RMF theories 
can realistically reproduce the bulk and single-particle properties of medium to heavy nu- 
clei. Successful RMF theories are characterized by large, neutral, Lorentz scalar and vector 
potentials (roughly several hundred MeV at nuclear equilibrium density), a new saturation 
mechanism for nuclear matter that arises from relativistic {i.e., velocity-dependent) inter- 
action effects from the scalar potential, and a nuclear spin-orbit force that is determined by 
this velocity dependence. The vector potential provides an efficient representation of the 
short-range nucleon-nucleon (NN) repulsion, while the scalar potential provides the mid- 
range NN attraction by simulating correlated, scalar-isoscalar two-pion exchange, which is 
the most important pionic contribution for describing the bulk properties of nuclear matter. 

For many years, RMF model parameters were determined by fitting them to the prop- 
erties of nuclei or nuclear matter using various "black-box" fitting schemes 0-0. Direct 
connections between the RMF parameters and the systematics of nuclear observables have 



begun to emerge only recently p|-p!T|. The modern theoretical viewpoint underlying the 
RMF, which relies on the ideas of Effective Field Theory (EFT) and of Density Functional 
Theory (DFT) [P!T|,P|,p!^ , provides a systematic truncation procedure for the RMF energy 
functional. However, it introduces more parameters than can be unambiguously determined 
by the relevant data. The present work uses EFT power counting to examine the quantity of 
information contained in the bulk and single-particle input data, and to identify how many 
parameters can be reliably determined by these data. 

The original Walecka model [|13| [without nonlinear field interactions), contained two 
free isoscalar parameters in nuclear matter, which were used to fit the equilibrium density 
and binding energy. An isovector coupling was added |I^ to reproduce the bulk symmetry 
energy, and for finite nuclei, a single scale parameter (the scalar "meson" mass) was chosen 
to reproduce the rms charge radius of "^^Ca. These procedures yielded reasonable results 
for charge-density distributions of spherical nuclei and for single-particle spectra, revealing 
that the nuclear shell model was obtained automatically without any specific adjustment 
of spin-orbit parameters [p!5|-[T7|,[5[| . Nevertheless, the nuclear compressibility and surface 



energy were too large in these simple calculations, leading to incorrect systematics in the 
total binding energies. Moreover, no quantitative measure of the accuracy of the results was 
obtained. 

It was also known at that time that the compressibility could be reduced by adding 
cubic and quartic scalar self-couplings to the model fl^. This procedure established a 
basic phenomenology for RMF studies of nuclei that is still used widely P^J5|J2(]|J^,pT|-pi 



Quantitative measures of the accuracy of the results were introduced, showing that there was 
a significant improvement over the original (linear) model. However, attempts to further 
improve the accuracy met with little success; for example, the introduction of a quartic 
vector self-coupling allowed one to perform RMF calculations of nuclei based on Dirac- 
Brueckner-Hartree-Fock calculations of infinite nuclear matter, but little improvement, if 
any, was achieved in reproducing the bulk and single-particle nuclear properties. 

An important advance was to define a measure of accuracy for the fit and then search the 
parameter space (in some particular model) to find a set of "optimal" parameters [|7|, ^ , |23 
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The conclusion of these surveys was that the hnear Walecka model plus cubic and quartic 
scalar self-couplings exhausted the capabilities of mean-field phenomenology to reproduce 
the relevant nuclear data 0]. However, the associated mean-field energy functional was trun- 
cated solely on the basis of its empirical success, essentially without theoretical justification. 
Our basic goal is to reassess this issue using the modern EFT framework for the nuclear 
many-body problem [pT] , ^ to ascertain how many model parameters (or linear combina- 
tions thereof) are accurately determined by input data drawn from bulk and single-particle 
nuclear properties. 



II. ENERGY FUNCTIONAL AND POWER COUNTING ESTIMATES 



Although we could use meson models as done historically, the analysis is more trans- 
parent with "point-coupling" models p7t^|r2| , which contain only nucleon fields in a local 
lagrangian. Because of the freedom to perform field redefinitions, a general point-coupling 
model is equivalent to a general meson model |jTl],|l2|. An energy functional of nucleon 



densities can be constructed by starting with a general point-coupling effective lagrangian, 
consistent with the underlying symmetries of QCD {e.g., Lorentz covariance, electromag- 
netic gauge invariance, and spontaneously broken chiral symmetry), and by constructing the 
corresponding one-loop energy functional. As discussed in Ref. this approach approx- 
imates a general DFT functional that incorporates many-body effects beyond the Hartree 
level when parameters are determined from finite-density data We expect this factor- 



ized representation of the functional to be reasonable because the large scalar and vector 
mean-field potentials imply that products of field operators in the lagrangian will be well 
approximated by products of densities in the functional ( "Hartree dominance" ) [|ll|,0 . More 



sophisticated approximations to the many-body problem {e.g., Dirac-Hartree-Fock, Dirac- 



Brueckner-Hartree-Fock, etc. [pO|-|35[|) could lead to more complete energy functional, but 
these calculations are not relevant for the questions studied here. 

The DFT framework generates a functional with an unlimited number of terms and cor- 
responding parameters. It is useful only if we can identify a valid expansion and truncation 
scheme. This requires an organization of terms in the effective lagrangian and a way to 
estimate the couplings. While precise relations between these couplings and the underly- 
ing QCD parameters are unknown, an estimate of the magnitude of the couplings can be 
obtained by applying Georgi and Manohar's Naive Dimensional Analysis (NDA) |]36,11,37 



The procedure is to extract from each term in the lagrangian the dependence on two 
primary physical scales of the effective theory: the pion decay constant ^ 94MeV and 
a larger mass scale A ^ ATrfj,/ ^JNj (where Nf is the number of light flavors) p8|j39[1 . The 
mass scale A is associated with the new physics beyond the pions: the non-Goldstone boson 
masses or the nucleon mass. This mass scale ranges from the scalar mass (~ 500 MeV) to 
the baryon mass (~ 1 GeV); the value A ^ 600 MeV is consistent with all existing empirical 
studies [11,12|- To establish the canonical normalization of the strongly interacting fields, 
an inverse factor of /tt is included for each field, and an overall factor of /^A^ fixes the 



normalization of the lagrangian. The physics of NDA is discussed further in Refs. ||Tl| and 



We restrict our attention to closed-shell nuclei, so that densities are static and three- 



vector currents are zero (see Ref. for generalized expressions). The energy functional is 
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therefore an expansion in powers of the nucleon scalar, vector, isovector-vector, tensor, and 
isovector-tensor densities scaled according to NDA: 



~ _ Ps _ V'^ 


(1) 


~ _ Pb _ V^V 

" /2A - /2A ' 


(2) 


^ P3 1 VVgV^ 

^'-/2A-2 /2A ' 


(3) 


_ Si ipa^'^i/j 

fA~ PA ' 


(4) 




(5) 



where if) denotes the nucleon field. The functional is also organized according to an expansion 
in powers of derivatives acting on these densities; NDA dictates that each gradient acting 
on a density is scaled by A: 

V^. (6) 

We will begin by considering only the isoscalar contributions to the energy functional, 
which dominate the bulk properties of nuclei. The functional is (following the notation of 
Ref. PI): 

/occ 
a 

+ fA^ J \^2ps - ^^diypsf + T^zPs + ^iPs + ^ipIPs + wip^ 

+ GPb - Cdiyp^f + CiPb - aiPs(Vps)^ - a2psiypBf 
+ isovector, tensor, and electromagnetic terms 

+ higher-order terms| . (7) 

Here the sum is over occupied single-particle states t/'q,, which may be different for protons 
and neutrons. Isovector and tensor terms will be discussed below. The electromagnetic 
energies follow closely the results from semi-empirical mass formulas . 



Naive dimensional analysis provides an organizational principle that directly translates 
into numerical estimates and an ordering of the terms in Eq. (J^). For example, each ad- 
ditional power of Ps is accompanied by a factor of /^A. The ratios of scalar and vector 
densities to this factor at nuclear matter equilibrium density are between 1/4 and 1/7 p7| , 
which serves as an expansion parameter. Similarly, one can anticipate good convergence for 
gradients of the densities, since the relevant scale for derivatives in finite nuclei should be 
roughly the nuclear surface thickness a, and so the dimensionless expansion parameter is 
1/Acr < 1/5. The expansion is useful because the coefficients have been shown empirically 
to be "natural," that is, of order unity PJl2 . 
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Since naturalness is valid, a semi-quantitative estimate of the contribution of each term 
in Eq. (|^) to the energy per particle can be made. In previous work, estimates were made 
for the contributions at the equilibrium point of symmetric nuclear matter by using the 
equilibrium density for both the baryon and scalar densities [|T0| , |12| . The estimates were 
compared to individual contributions directly evaluated from Eq. (^, but many terms do 
not contribute in nuclear matter. 

Here we extend this analysis to finite nuclei with Z protons, iV neutrons, and A = A^ + Z. 
This requires estimates for the isoscalar scalar and vector densities in the finite system as 
well as for isovector, tensor, and gradient densities. Since we must allow for variations 
in natural coefficients as well as uncertainties in the scale A, it is sufficient to use simple 
local-density approximations. 

The scalar and vector densities differ by less than 10% in ordinary nuclei, so we can use 
the average baryon density for both: 

Ps,Pb — ' TiT ' (8) 

where the average value of /(x) is calculated as 

l^f{x))^^^d'xi{x)p^{x) . (9) 

Like the surface thickness, the average nuclear density (pe) is essentially model independent. 
We determine it for each nucleus by using densities from one of the best-fit models (all give 
the same results to within a few percent). 

A local-density estimate (per nucleon) for a general isoscalar term in Eq. (|^ with di- 
mensionless coefficient /3 is made according to 

jfA'jd'xpip^np^y ^ PA (^^y ' , (10) 

and we take l/2</5<2asa reasonable range of natural coefficients. These estimates 
with the associated error bars from the variation in (3 are shown as small squares in Fig. [1| 
for ^^O and in Fig. ^ for ^^'^Pb. The magnitudes of energy contributions in Eq. (|^) from 
two representative RMF point-coupling models from Ref. are shown as larger unfilled 
symbols (one model on each side of the error bars). These models provide very accurate 
predictions of bulk nuclear properties |T^. The energy contributions are determined for each 



nucleus by making multiple runs while varying each parameter slightly around its optimized 
value, which enables us to deduce the logarithmic derivative with respect to each parameter. 
The filled symbols denote the sum of the values for each power of the density. The binding 
energy per nucleon in nuclear matter is denoted with eo. 

The representative RMF models validate the isoscalar estimates from Eq. (|l3), and the 
resulting hierarchy of isoscalar contributions is quite clear. The question is: How far down 
in the hierarchy can we reliably determine contributions and their associated parameters? 
To address this question, we can consider a series of models truncated at different levels in 
the hierarchy. As additional terms are added, the reproduction of experimental observables 
should improve. When improvement stops, contributions beyond this point are unresolved 
and irrelevant. 
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FIG. 1. Contributions to the energy per particle in ^^O determined by logarithmic derivatives 
with respect to the model parameters (see text) for two RMF point-coupling models from Ref. jl2| . 
Absolute values are shown. The filled symbols are net values. The small symbols indicate estimates 
based on NDA [Eq. (0)], with the error bars corresponding to natural coefficients from 1/2 to 2. 
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FIG. 2. Same as Fig. | for ^ospb. 
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FIG. 3. Figure of merit ("x^") values for RMF meson models [^:12], ordered according to the 
truncation level of the energy functional. 



To make a meaningful analysis, we must decide on the most relevant observables. The 
Density Functional Theory (DFT) perspective motivates the observables that we expect to 
calculate reliably with a RMF model [|T|. These include the total binding energy, density 
distributions at low momentum transfer (including rms radii), the chemical potential, and 
splittings between quasi-particle levels near the Fermi surface. We construct a measure of 
the accuracy of the fit to these observables, which is based on squared deviations between 
predictions and data for various observables but is not a statistical x^- Instead, we apply 
relative weightings based on our DFT expectations, with the total energy weighted most 
heavily. Properties of five doubly magic nuclei were used (see Refs. and [l^ for details). 

In Fig. 1^, the impact of different RMF meson model truncations is shown by plotting this 
figure of merit ( "x^" ) against the maximum power of fields in a truncated energy functional. 
We have performed this test with both point-coupling and meson models, with similar 
conclusions. We focus here on the meson model results because we are able to directly 
compare them to standard phenomenological models. (See Ref. [O] for the point- coupling 
analog.) The "full" models (which include all nonredundant terms at a given order) show 
that one needs to go to the fourth power of fields to get the best fits, but going further yields 
no improvement. Analogous behavior is found for RMF point-coupling models with powers 
of densities replacing powers of fields, although the quality of the fit is much better when 
the truncation is at the cubic level 
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Thus contributions to the energy/particle at the 
level of 1 MeV or so are at the limit of resolution. Fifth-order isoscalar contributions to the 
energy/particle, which are predicted to be less than 1 MeV, are simply not determined by 



the optimization [|T2 

The "0" only" results in Fig. ^, which include only scalar fields for n > 2, show that nearly 
optimal results can be obtained with just a subset of terms at each order. This explains 
the successes of the most widely used RMF models, which add only 0^ and 0^ terms to the 
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FIG. 4. Coefficients for four accurately fit RMF point-coupling models from Ref. Each 
model is represented by a different shape, and the shading shows the type of term (scalar, vector, 
or mixed). 
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TABLE I. Improved coefficients for point-coupling RMF models [^] . 



coefficient linear combination density scaling deduced value 

0,1 K2 + C2 p+ -0.51 ±0.01 

^3 K3 + m pI +1.3 ±0.1 

^^2 ^2jC2 P- -2.0 ±0.4 

K4 + Ci + m pI -2.4 ±0.7 

K3-^?i/3 p+p- ±0.2 ±1.0 

^6 K4-C4 pIp- -2.6 ±0.8 

K4 ± C4 - ??22/3 p+p'i 



original Walecka model lagrangian. Since the coefficients are underdetermined by the data, 
it is not surprising that the systematic optimization of more complete RMF models yields 
multiple EFT parameter sets with accurate reproductions of nuclear properties [11,12||- '^^^ 



variation of coefficient values in Eq. (|^ provides a measure of how well the parameters are 
actually determined by the data. In Fig. ^, the seven coefficients of isoscalar, non-gradient 
terms from four point-coupling models are plotted. The analogous results for coefficients 
in several meson models is shown in Fig. |^. We note that all coefficients are natural, i.e., 
order unity. However, the spread in coefficient values is significant, particularly for the 
meson models, and does not correspond well to the power-counting order. We conclude 
that different linear combinations of the coefficients must be considered to draw reliable 
conclusions about how many are determined by the data. 

Can we find a more systematic power counting scheme? The similar size of the scalar 
density ps and the vector density pe suggests that we count instead powers of p+ = (ps±Pb)/2 
and p_ = (ps — Pb)/2. The corresponding "improved" coefficients are listed in Table | [see 



Eq. (13)]. The spread in these coefficients for four RMF point-coupling models from Ref. ||T2 
are shown in Fig. ^ The terms are organized according to the power of p+ and p_, with p_ 
behaving as p^^ [T^. 

The leading orders are very well determined, with a systematic increase in uncertainty. 
Even the sign is undetermined for the parameter ^^4, which is shown with unfilled symbols, 
but the next parameter (Qq) appears to be reasonably well determined. Higher-order terms 
are undetermined from the optimizations. Deduced values and uncertainties based on this 
sample of models are given in Table |. We see that of the seven isoscalar, non-gradient 
parameters in Eq. (^, four linear combinations are clearly determined by bulk nuclear 
observables, with probably a fifth combination as well. 

Now we return to consider a more complete version of Eq. (0) by including the leading 
isovector and tensor terms: 

/occ 
a 

+ fl^^ J {^2ps - fidi^psf + Ksp^ + Kip^ ± r/iPaPs ± f]2plps 

± GPb - Ci(VpB)' ± C4Pb - 5lPs(Vps)' - 52Ps(VPb)' 

+ i2pi - ld(Vp3)^ ± tippIp^ ± /„Vpb ■ s ± /pVps ■ 
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FIG. 6. Improved coefficients for tlie same four models as in Fig. The "order" is determined 
by counting powers of p+ and p_. 

+ electromagnetic terms + higher order termsj . (11) 

We will need estimates for isovector, gradient, and tensor densities. An estimate for the 
isovector density is obtained by a simple scaling of the local-density approximation to the 
baryon density: 



P3 



flk 



(12) 



which reproduces the direct average of isovector densities in best-fit models to about 10% 
accuracy. We use [see Eq. 



(13) 



for gradients and 



(14) 



for tensor terms. 

Estimates for all of the terms in Eq. ([TT|) and the corresponding contributions from the 
two point-coupling models are indicated with different small symbols in Figs. |^ (^^O) and 
g ^208p^^_ Solid inverted triangles represent the net contribution of two gradient terms at 
second and third order in the density expansion. The actual contributions from the direct fits 
for the two nuclei from isovector, gradient, and tensor terms exhibit significant differences, 
but the estimates follow these systematics quite reasonably. This enables us to predict how 
many parameters of each type can be determined reliably. 

The isovector terms appear only on the graphs for ^°^Pb. The factor (A^ — Z)/2A, which 
is only 10% even for Pb, severely cuts down the sensitivity to isovector terms (especially 
since the factor must appear with even powers). The magnitude of the leading isovector 
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term (oc ^2pi ) is comparable to the fourth-order isoscalar term, which is at the hmit of 
what can be determined rehably from fitting the binding energy. Fits of the subleading 
term (oc fjppsp^) are border hne unnatural, which is typical for energy contributions that 
are essentially unresolved. Furthermore, while individual values for ^2 and rjp vary widely 
between different parameter sets, the combination ^2 + VpiPs), with (ps) ~ (Pb) ~ 0.16 
from averaging over nuclei from oxygen to lead, is well determined. We conclude that only 
one isovector parameter is determined by the bulk observables. A term proportional to 
(ps)^, which one might think important for the neutron-star equation of state, is completely 
undetermined by ordinary nuclei.Q Our results support previous observations in meson 
models that at the mean-field level, the role of the p meson is simply to reproduce the 



correct symmetry energy Including any other isovector degrees of freedom, such as 

an isovector, scalar meson, is unnecessary. However, the use of specific isovector observables, 
such as the difference between neutron and proton distributions at low momentum transfer, 
may provide more constraints on the isovector parameters. 

The energy estimates for isoscalar tensor terms imply that only one parameter, at best, 
can be determined. Higher-spin terms, which will require more gradients and have smaller 
average densities, are not at all constrained. The tensor terms are interesting to consider 
because a fraction of the spin-orbit force can be generated by including an isoscalar tensor 
coupling of the vector field to the nucleon. At the same time, the scalar and vector fields will 
be reduced (and therefore the effective nucleon mass M* will be increased). Nevertheless, the 
spin-orbit potential arises predominantly from the large scalar and vector fields; attributing 
more than one-third of the potential to the tensor coupling produces unrealistic surface 



systematics In Ref. |^6|, it was shown how to accurately relate the size of the tensor 
coupling to the value of M* in nuclear matter and the size of spin-orbit splittings in finite 
nuclei. 

Finally, the gradient terms follow the same pattern in the energy: the leading term 
is barely above the limit of resolution. In fact, there are two isoscalar gradient terms at 
leading order (scalar and vector), but if we consider linear combinations based on p+ and 
p_ (see the improved coefficients Aj and Ai in [1^), only their sum is well determined. The 
subleading-order contributions appear to be highly unnatural individually, but their sum 
almost vanishes. Thus we conclude that only one gradient parameter is determined. This 
has implications for the possibility of resolving effective meson masses from fits of point- 
coupling parameters. In particular, estimates for the scalar and vector masses squared are 
A^K2/Kd and h?C,2/Cd, respectively. Since Kd and Q are not individually determined, only 
the difference of the inverse squared masses can be extracted. 



III. IMPLICATIONS AND SUMMARY 

We have shown that EFT power counting and naturalness lead to reliable estimates 
of the energy contribution in a finite nucleus from any term in a general RMF energy 



Fortunately, the condition of beta equilibrium in a neutron star greatly restricts the role of a 
(^3)^ term @. 
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functional. By examining energy systematics based on these estimates and the spread among 
coefficient values in models fit to a given set of nuclear observables, we can deduce how many 
parameters are actually constrained by those data. We conclude that four (possibly five) 
isoscalar, non-gradient parameters, one gradient parameter, and one isovector parameter 
are well determined by the usual bulk nuclear observables (binding energies, charge density 
distributions, and spin-orbit splittings in doubly magic nuclei). Since conventional RMF 
models include this distribution of parameters (plus two additional unnecessary gradient 
terms), it is not surprising that empirical surveys have found no phenomenological reason 
to add additional terms 0. 

The handful of parameters that are well determined can be associated with an equal 
number of nuclear properties and general features of RMF models. In particular, 

1. Two isoscalar, non-gradient parameters are very well determined. These correspond 
to the highly constrained values for the equilibrium density {kp ~ 1.30 ± 0.01 fm^ ) 
and binding energy (16.0 ±0.1 MeV) of observed nuclear matter. 

2. An additional isoscalar constraint is that M* ^ 0.61 ±0.03, if the isoscalar tensor term 
is set to zero. This range ensures an accurate reproduction of spin-orbit splittings in 
finite nuclei. Small increases in M* without changing the splittings can be accom- 
modated by including an isoscalar tensor term, with the trade-off given by a simple 
local-density approximation in Ref. [HO . 



3. A fourth isoscalar constraint comes from the nuclear matter compressibility. The 
constraint is much weaker, in the range of ~ 250 ± 50 MeV. 



4. The possibility of a fifth isoscalar constraint has been considered by Gmuca |2^, who 
argued that one needed separate scalar and vector fourth-order terms to tune the 
density dependence of the scalar and vector parts of the baryon self-energies. This 
would correspond to constraining Qq from Table |. In addition, some form of isoscalar 
nonlinear vector coupling is needed to soften the high-density equation of state to be 
consistent with observed neutron star masses 



Since only one isoscalar gradient parameter is determined, it is not useful to allow 
the scalar and vector masses (or their equivalents in a point-coupling model) to vary 
independently. Thus it is convenient to fix the vector mass at a natural size, such as 
the experimental mass for the u. Then a scalar mass of 500 ± 20 MeV is required. (It 
may be fine tuned using the ^°Ca rms charge radius 3.45 fm as a target.) 

The one isovector parameter can be fixed by the surface-corrected volume symmetry 



energy which falls in the range 34 ± 4 MeV 0. Since no isovector gradient is 
determined, setting the isovector vector meson mass to the experimental p meson 
mass is adequate. 

To produce the best quantitative fits to nuclear properties for any given model, the model 
parameters should be obtained from an optimization procedure involving calculations of a set 
of finite nuclei from oxygen to lead |T0| , p!2| . However, to ensure a reasonable (if not optimal) 



reproduction of finite nuclei, it is sufficient to reproduce the nuclear matter properties given 
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above. Note that all properties must be satisfied and the resulting parameters must be 
natural.3 

This last point has important implications for an alternative class of RMF investiga- 
tions, which are not based on EFT but are defined by some specific underlying physics. 
Examples include models that contain different degrees of freedom (quarks) p7| , ^ , or dif- 
ferent types of meson-baryon couplings (derivative couphngs) p9| , ^ , or explicit inclusion 
of chiral symmetry |pT|-|5B[]. (The chronology for this development can be gleaned from some 
recent review articles the papers cited above, and references therein.) Unfortunately, 

in many of these studies it was believed to be sufficient to reproduce just the nuclear matter 
equilibrium point, and thus the results for finite nuclei were not as accurate as expected 
from the EFT analysis |]9|,|46||. One cannot justify the underlying physics of a model if it can 
reproduce only a subset of the nuclear calibration data. 

Moreover, both the input observables and the parameters are highly correlated and noisy; 
therefore, it is inappropriate to use the deduced parameters from Table |I| instead of a direct 
fit to observables. It is possible to achieve a sensible picture of nuclei by this procedure but 
not a very accurate fit. 

The improved coefficients from Table | suggest that a nonrelativistic point-coupling EFT, 
with an expansion in only, should be consistent with similar power-counting estimates. 
Indeed, the phenomenologically successful Skyrme models are in this category, and an NDA 
study shows that they are natural ||5^. A scale of A = 600 MeV is consistent with the trends 
in all of the relativistic and nonrelativistic models. A comparison of these approaches from 
the EFT perspective will be given elsewhere. 

The EFT/DFT perspective provides insight into which additional observables might fur- 
ther constrain the parameters of a general RMF functional. Since four or five isoscalar, 
non-gradient terms in the energy functional are determined, it is unlikely that additional 
bulk, isoscalar observables will provide further constraints. This is consistent with calcu- 
lations of deformed even-even nuclei, which reproduce experimental systematics once the 
energy functional is constrained to reproduce bulk properties of doubly magic nuclei. On 
the other hand, since only one isovector parameter is fixed, additional isovector observables 
could be relevant. In particular, an accurate measurement of the neutron radius in a heavy 
nucleus, which would result from proposed parity violation experiments, may fix an ad- 
ditional isovector parameter. Additional constraints on gradient terms might be provided 
from data on ground-state currents at nonzero (but low) momentum transfer, for example, 
from elastic magnetic scattering. An assessment of these constraints within the framework 
described here is in progress. 
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further caution is that there are many correlations among these properties, so that the allowed 
ranges should not be considered truly independent. 
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